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The  plasma  plume  from  a  200  W  Hall  Effect  Thruster  fueled  by  iodine  vapor  was 
analyzed.  The  plasma  source  included  a  laboratory  propellant  feed  system  and  a  laboratory 
model  Hall  thruster  powered  by  a  breadboard  power  processing  unit.  The  hollow  cathode 
was  fed  with  xenon.  The  distribution  of  iodine  ions  was  measured  with  an  ExB  probe,  an 
electrostatic  analyzer  (ESA),  and  a  combined  ESA/ExB  probe.  The  distribution  of  xenon 
ions  was  also  measured.  Multiply  charge  species  were  detected  with  both  iodine  and  xenon. 
Significant  populations  of  diatomic  iodine  ions  were  also  detected.  The  dimer  fraction  was 
found  to  vary  with  operating  conditions  and  angular  distance  from  the  beam  centroid.  The 
greatest  populations  of  high  energy  dimers  were  observed  at  lower  discharge  current,  and  at 
angles  of  30-50  degrees  from  the  beam  centroid.  At  one  operating  condition,  the  high  energy 
dimer  fraction  was  approximately  20%  by  mass  in  this  region.  Significant  dimer  fractions 
could  increase  thrust  to  power  by  useful  amounts  at  the  cost  of  some  overall  efficiency. 
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Nomenclature 

magnetic  field 
electric  field 

charge  of  an  electron,  1.6  x  10"19  C 
beam  momentum  in  axial  direction 
species  fraction 

gravitational  constant  at  Earth’s  surface,  9.81  m/s2 

current,  subscripts  b  for  beam,  c  for  cathode,  d  for  discharge,  m  for  magnet 
specific  impulse 

current  density 

mass  flow  rate  divided  by  discharge  current 
characteristic  distance/dimension 

mass  flow  rate,  subscripts  a  for  anode,  c  for  cathode,  i  for  beam  ions 
ion  mass 
discharge  power 

pressure,  subscript  s  for  sensor,  x  for  xenon 
ion  charge 

radial  direction  or  distance 

thrust 

time 

discharge  potential 
velocity 

mean  ion  charge 

axial  direction  or  distance 

efficiency 

flux 

viscosity 


1  Chief  Scientist  for  Hall  Thrusters,  AIAA  Associate  Fellow 

2  Research  Engineer,  AIAA  Member 

1 

American  Institute  of  Aeronautics  and  Astronautics 
Approved  for  public  release;  distribution  unlimited. 


I.  Introduction 


THE  Hall  Effect  Thruster  (HET)  is  an  efficient  form  of  electric  propulsion  that  is  typically  used  in  space  for 
orbit-raising  and  orbit  maintenance.  Hall  thrusters  were  invented  in  the  United  States,  but  first  transitioned  to 
flight  in  the  Soviet  Union.  1  2'  ’  4  The  first  American  HET  to  fly  in  space  was  the  Busek  BHT-200,  launched  in 
2006  as  part  of  the  US  Air  Force  TacSat-2  satellite  and  fueled  by  xenon.5  Testing  of  a  BHT-200  fueled  by  iodine 
vapor  yielded  stability  and  performance  comparable  to  that  observed  with  xenon.6  Plume  current  measurements 
with  the  same  thruster  showed  lower  beam  divergence  and  the  presence  of  E+  near  the  beam  centroid.7  Testing  of  a 
1-kW  thruster  fueled  by  iodine  also  showed  high  efficiency  and  low  beam  divergence.8 

This  paper  report  iodine  ion  energy  and  species  distributions  measurements  carried  out  with  an  experimental 
version  of  the  BHT-200  shown  in  .  These  measurements  were  obtained  with  an  ExB  probe,  an  electrostatic  analyzer 
(ESA),  and  a  combined  ESA/ExB  probe  at  various  operating  conditions  and  angles  with  respect  to  the  beam 
centroid.  The  paper  also  provides  a  discussion  of  how  the  composition  of  the  iodine  plume  may  affect  efficiency 
and  thrust  to  power. 

A.  Hall  Thruster  Physics 

A  Hall  effect  thruster  uses  crossed  electric  and  magnetic 
fields  to  generate  and  accelerate  ions.  The  gas  in  the  discharge  is 
partially  ionized,  although  propellant  utilization  may  approach 
100%.  The  overall  structure  of  the  BHT-200  is  defined  by  a 
magnetic  circuit  that  produces  a  steady  magnetic  field,  B  ,  in  the 
nominal  radial  direction  across  an  annular  channel.'’  Neutral 
propellant  is  introduced  through  the  anode  at  the  base  of  the 
channel.  The  downstream  portion  of  the  channel  is  formed  by  a 
dielectric  material,  between  which  the  bulk  of  the  plasma 
discharge  occurs.  A  potential  difference,  Vd ,  is  applied  between 
the  anode  and  a  hollow  cathode  located  outside  the  channel.  The 
resulting  electric  field,  E ,  is  predominantly  axial  and  is 
concentrated  near  the  channel  exit  by  interactions  between  the 
magnetic  field  and  the  plasma.  Typically  the  cathode  is  allowed 
to  float  with  respect  to  facility  ground. 

In  the  channel,  electrons  are  strongly  magnetized  and  their 
transport  is  predominantly  azimuthal  due  to  the  ExB  Hall  Effect.  The  extended  electron  path  enables  an  efficient, 
impact  driven  ionization  cascade.  The  peak  electron  temperature  is  typically  several  tens  of  volts.  Axial  transport 

across  B  is  provided  by  collisions,  wall  effects,  and  additional,  “anomalous”  transport  identified  in  the  early 
1960s.10  Ions  are  weakly  magnetized  and  most  are  accelerated  directly  out  of  the  channel,  forming  the  ion  beam. 
The  bulk  of  the  cathode  electrons  go  to  neutralize  the  ion  beam  current,  I h  .  However,  some  circulate  back  through 
the  channel,  driven  by  the  potential  difference.  These  electrons  seed  the  discharge  cascade.  A  control  volume 
analysis  shows  that  the  discharge  current,  I d  ,  is  equal  to  the  cathode  current,  Ic ,  which  is  equal  to  the  sum  of  I  h 
and  the  electron  flux  to  the  discharge. 

B.  Rationale  for  Iodine 

Iodine  propellant  may  be  used  in  Hall  thrusters,  ion  engines,  and  other  types  of  plasma  based  thrusters.  Table  1 
compares  some  key  properties  of  iodine  with  those  of  the  noble  gas  xenon.  Atomic  iodine  is  lighter  than  xenon,  but 
iodine  stores  and  vaporizes  as  a  diatomic  molecule,  I2.  Furthermore,  ionization  energies  and  electron  impact 
ionization  cross  sections  for  I  and  I2  are  larger  than  ionization  energies  and  cross  sections  for  Xe.1411  These 
attributes  make  iodine  easier  to  ionize  than  xenon,  which  result  in  similar  or  greater  propellant  utilization.  Iodine 
performance  measurements  indicate  little  or  no  efficiency  penalty  with  respect  to  xenon.6  8  At  equivalent  operating 
points,  thrust  to  power  sometimes  increases,  but  this  is  typically  counterbalanced  by  lower  specific  impulse. 
Measurements  also  show  consistently  lower  beam  divergence  with  iodine. 


Fig.  1  Iodine  Thruster  in  Test  Facility 
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Table  1  Properties  of  iodine  and  xenon.12  13  14 


Element 

1 

Xe 

Atomic  Mass 

126.9 

131.3 

Ionization  Properties  (monatomic) 

First  Ionization  Potential  (eV) 

10.5 

12.1 

Peak  Cross  Section  (10"16  cm2) 

6.0 

4.8 

Storage  and  Handling  Properties 

Storage  density  (gm/cm3)  near  room  temp. 

4.9 

1.6* 

Melting  Point  (°C) 

113.7 

-112 

Boiling  Point  at  10  Pa  (°C) 

9 

-181 

*14  Mpa,  50  C  (NIST  Database) 

Storage  and  handling  properties  for  iodine  and  xenon  are  very  different.  Iodine  stores  in  the  solid  phase  at 
approximately  three  times  the  density  of  xenon.  The  vapor  pressure  required  to  fuel  a  thruster  (much  less  than  1 
atmosphere),  may  be  generated  by  heating  the  solid  to  a  modest  temperature,  e.g.  80  -  100  degrees  C.  Xenon  must 
be  stored  in  h  high  pressure  tanks  or  at  cryogenic  conditions. 

The  only  significant  disadvantage  to  iodine  is  that  it  is  a  halogen,  which  means  it  is  reactive  with  some 
substances.  However,  the  vapor  pressure  of  U  is  greater  than  1  Pa  at  -12  C,  which  suggests  minimal  deposition  on 
spacecraft  surfaces.  Furthermore,  inside  the  thruster,  condensation  will  not  occur  at  typical  operating  temperatures. 
Reactivity  with  wetted  surfaces  is  an  issue,  but  this  may  be  addressed  though  material  selection. 

C.  Unique  Features  of  the  Iodine  Discharge 

In  a  typical  Hall  thruster,  ionization  is  produced  through  inelastic  scattering  of  high  energy  electrons  off  a  neutral 
gas  that  is  flowing  through  the  channel.  The  electron  impact  ionization  reactions  for  iodine  are  as  follows: 


e  +  7  — >■  2e  +  /+  (la) 

e  +  I  — »3e  +  /2+  (lb) 

c  + I2  — ^  2t?  +  /->  (lc) 

e  +  I2  — »  3e  +  I2+  (Id) 


Fig.  2  shows  electron  impact  ionization  cross  sections  for  these  reactions.15  The  ionization  energies  for  I  and  I2  are 
10.5  and  9.3  eV,  respectively. 


Electron  Energy  (eV) 


Fig.  2  Iodine  ionization  cross  sections 
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Dissociation  is  also  caused  by  inelastic  scattering: 

£?  +  /■>  — ^  {?  +  /  +  /.  (2a) 

Although  the  I2  dissociation  energy  of  1.57  eV12  is  far  lower  than  the  electron  temperature  in  the  discharge,  the 
measured  presence  of  in  the  plume  suggests  dissociation  is  incomplete.  Dissociation  processes  may  also  produce 
ions  through  the  reactions 

c  + 1 2  — ^  +  /  + 1+  (2b) 

and 

e  +  I2  — »  3e  +  /+  +/+  .  (2c) 

Dissociative  ionization  obviously  requires  more  energy  than  mere  dissociation,  but  electron  temperatures  of  ~20  eV 
are  predicted  in  the  channel  of  a  xenon  fueled  BHT-200  by  the  hybrid  Particle-in-Cell  modeling  code  HPHall. 1617 

Negative  ions  may  also  be  present  in  the  iodine  discharge.  These  are  formed  through  attachment  processes.  If 
formed,  negative  ions  would  move  upstream  toward  the  anode,  where  they  would  most  likely  be  recycled 

Charge  exchange  collisions  are  also  present.  These  are  particularly  important  in  the  ion  beam,  where  the 
particles  have  long  distances  to  interact.  Some  reactions  include 


I  + 1+  ~^I+  +/, 

(3a) 

/+/2+^r+/2, 

(3b) 

/+/2+  h>/+  +/+, 

(3c) 

i2+i+  ^/2++/, 

(3d) 

/2  +  4+  — » I2  +  4, 

(3e) 

1 2  +I2+  ->/2  +/+. 

(3f) 

All  of  the  listed  ionic  species  are  measured  in  the  plume. 

Over  long  distances,  ion  recombination  could  significantly  diminish  the  dimer  population  and  affect 
interpretation  of  plume  measurements.  However,  recombination  is  likely  to  be  minimal,  and  in  modeling  xenon  Hall 
thrusters,18  it  is  often  ignored. 

In  prior  work,  the  composition  of  the  flight  model  BHT-200  thruster  plume  was  measured  with  an  ExB  probe 
located  near  the  beam  centroid  at  one  operating  condition  ( Vd  =  250  V,  Pd  -  200  W).  The  primary  iodine  ionic 
species  in  the  plume  are  I+,  I  ,  and  I2+.  Per  Table  2,  approximately  3%  of  the  iodine  plume  was  found  to  consist  of 
diatomic  ions.  This  was  twice  the  size  of  the  measured  l21  fraction  and  40%  larger  than  the  Xe2+  fraction.6'7 
However,  variations  in  space  and  with  operating  conditions  were  not  explored.  The  degree  of  dissociation  may  not 
be  uniform  throughout  the  discharge  and  plume,  or  constant  with  respect  to  plasma  density  and  electron  temperature, 
which  vary  with  Id  and  Vd  . 

Table  2  Iodine  molar  flux  fractions  at  250  V  and  -200  W 


Species 

Mole  Fraction  (%) 

l+ 

95.3 

l2+ 

2.9 

l2+ 

1.5 

r 

0.3 

D.  Kinetic  View  of  the  Beam 


The  beam  ion  current  is  the  sum  of  k  ionic  species: 


4  = 


X7*  =Xw* 

k  k 


9  k 
Mk 


(4) 


Here,  q  is  the  ion  charge,  M  is  the  ion  mass,  and  m  is  the  total  species  mass  flow  rate.  Thrust  is  determined  by  the 
axial  (Z-direction)  momentum  carried  by  the  heavy  particles 
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(5) 


F  =  X  mk  <  Vk,z  >• 

k 

Most  momentum  is  carried  by  the  ionic  species,  although  modeling  and  testing  of  xenon  thrusters  shows  a  small  but 
measurable  fraction  of  the  momentum  is  carried  by  neutrals.  The  momentum  carried  by  an  individual  ion  at  energy 
e  is  given  by  p  =  Mv ,  where 


v  = 


(6) 


Here,  V  is  the  potential  through  which  the  ion  has  been  accelerated,  given  by  V  =  Vd  -  /  A  Vi  .  The  loss  term 
accounts  for  the  floating  cathode  potential  and  the  point  of  ion  origin  with  respect  to  the  anode  potential. 

The  mean  exhaust  velocity  or  specific  impulse  is  defined  by  thrust,  Isp=F/mg0.  Here,  m  is  the  total  thruster 
mass  flow  rate,  and  g0  is  the  force  of  gravity  at  the  surface  of  the  Earth.  The  ratio  of  thrust  to  power  (F  /  P) 
available  at  a  particular  I  determines  the  efficiency: 


F  Iygp  F2 
P  2  2  mP 


The  relationship  between  FtP  and  /  is  plotted  for  thruster  efficiencies  between  50%  and  100%.  The  only  way 
to  increase  F  /  P  at  constant  /  is  to  increase  efficiency. 


Specific  Impulse  (s) 


- 50% 

—  •  -60% 
- 70% 

----80% 

. 90% 

100% 


Fig.  3  Relationship  Between  Thrust  to  Power  and  Specific  Impulse 


For  thruster  analysis,  often  only  the  anode  propellant  flow  rate  and  discharge  power  are  considered  when 
formulating  the  efficiency.  This  anode  or  thrust  efficiency  may  be  expressed  as  the  product  of  the  propellant 
utilization  efficiency  ( iju ),  electrical  efficiency  ( r\e ),  acceleration  efficiency  (rja)  and  beam  efficiency,  (r/b ), 

Tlt=TluT1eT1aTlb-  (8) 

The  propellant  utilization  efficiency  is  the  ratio  of  the  beam  ion  mass  flow  rate  to  the  anode  mass  flow  rate; 
iju  =  rhf  /  ma  .  The  electrical  efficiency  is  the  ratio  of  ion  beam  current  to  discharge  current,  r]e  =  Ibl  I d  .  The 
acceleration  efficiency  is  the  ratio  of  the  total  potential  through  which  the  ions  are  accelerated  to  the  discharge 
potential;  r]a  =(V IVd)  .  Following  the  method  of  Ref.  19,  the  beam  efficiency  may  be  expressed  as  rib=rjsFt2. 
The  thrust  correction  factor,  Ft ,  accounts  for  momentum  lost  to  beam  divergence.  The  species  efficiency,  rjs , 

5 

American  Institute  of  Aeronautics  and  Astronautics 
Approved  for  public  release;  distribution  unlimited. 


accounts  for  the  momentum  spread  in  a  beam  with  particles  of  multiple  sizes.  When  analyzing  colloidal  thrusters, 
Lozano  called  this  concept  polydispersivity. 20  This  concept  was  referred  to  as  propellant  efficiency  by  Brown  et 
al.21 


If  the  thruster  efflux  is  assumed  to  be  symmetric  about  the  beam  centroid,  the  number  of  particles  crossing  a 
hemispherical  surface  distance  r  from  the  thruster  exit  is  given  by 


dN  t  f 

- =  2 w2  I  V(9)sm(9)d9. 

dt  J 


W2 


(9) 


Here,  9  is  measured  from  the  plume  centroid,  and  T(9)=n<vr  >  is  the  local  flux  in  the  radial  direction. 
Integrating  the  ion  current  density,  j(9) ,  over  the  surface  gives  the  beam  current; 


nil 


Ih=2m-2  ^  j{9)s'm(9)d9  .  (10) 

o 

By  species,  the  molar  (number)  fraction  of  r(<9)  crossing  the  surface  is  given  by  hk ,  where  ^hk  =  1 .  Thus,  ion 

k 

current  by  species  is  proportional  to  ^  hkqk  and  mass  flux  by  species  is  proportional  to  ^  hkMk  .  Here,  the 

k  k 

normalized  current  and  mass  flux  fractions  are  called  fk  and  gk ,  respectively. 


The  thrust,  which  is  equivalent  to  the  change  in  axial  momentum  per  unit  time,  is  given  by 


71 1 2 


F  =  — =  2/zr2  f  TMvr  cos{6)sm{0)d6  . 
dt  J 


(ID 


If  the  momentum  flux,  TMvr ,  is  assumed  to  vary  in  proportion  to  j(9) ,  then  the  fraction  of  the  beam  momentum 
that  is  in  the  axial  direction  may  be  estimated  through  the  relation 


7T  I  2 

2m-2  J /((9)sin((9)cos(^)a,(9 

Ft= - 5 - - - •  (12) 

h 

This  relationship  essentially  assumes  that  average  ion  speed,  charge,  and  mass  are  constant  with  angle.  However,  as 
reported  in  this  paper,  hk  ,  fk  and  gk  may  vary  across  the  plume.  Furthermore,  ion  energy  measurements  here  and 


historically  show  non-Maxwellian  ion  velocity 
distributions  that  vary  with  proximity  to  the  beam 
centroid  and  thruster  exit.2* 

The  degree  to  which  the  local  composition  of 
the  plume  affects  thruster  performance  may  still 
be  gauged  by  j{9) ,  which  decreases  with  9  at  a 
roughly  exponential  rate.  However,  on  doing  so, 
the  weighting  factors  of  equations  9-12  must  be 
considered.  The  flux  and  beam  current  weighting 
factor,  sinff?) ,  peaks  at  90  degrees,  while  the 
thrust  weighting  factor,  sinff?)  cos(f?) ,  peaks  at  45 
degrees.  This  counteracts,  to  some  extent,  the 
exponential  decay  of  the  plume  and  it  makes  the 
composition  of  the  beam  out  to  at  least  45  degrees 
significant  with  regard  to  performance.  The 
weighting  functions  are  plotted  in  Fig.  4. 


Fig.  4  Weighting  factors  as  a  function  of  angle 
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II.  Apparatus  and  Procedure 

In  this  research,  the  plume  of  a  low  power  laboratory  model  iodine  fueled  Hall  thruster  was  measured  with 
electrostatic  and  electromagnetic  plasma  wetted  probes. 

A.  Hall  Thruster  System 

The  Hall  thruster  used  for  this  research  was  an  experimental  version  of  the  Busek  BHT-200.  The  beam  was 
neutralized  by  a  barium  impregnated  hollow  cathode  flowing  xenon.  The  system  was  powered  by  a  breadboard 
power  processing  unit  (PPU).  Operationally,  the  thruster  was  pre-heated  with  a  xenon  discharge  before  switching 
over  to  iodine.  Fig.  5  shows  the  iodine  plume  during  testing. 


B.  Xenon  Feed  System 

Xenon  flow  was  regulated  by  Unit  7300  series  flow  controllers  sourced  from  Unit  Instruments.  The  anode  was 
fed  with  a  50  SCCM  controller,  which  was  calibrated  to  minimize  uncertainty  at  flow  rates  of  interest.  The  cathode 
was  fed  with  a  10  SCCM  controller. 

C.  Iodine  Feed  System 

I2  was  sublimed  from  high  purity  crystals  inside  a  thermally  controlled  reservoir  located  outside  the  vacuum 
tank.  The  I2  vapor  passed  through  a  heated  line  to  the  thruster  located  under  vacuum.  The  flow  rate  was  controlled 
by  regulating  the  temperature  (hence  pressure)  at  the  source  as  described  in  Ref.  6.  The  system  is  shown  in  Fig.  6. 
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D.  Test  Facility 


The  system  was  tested  in  Busek’s  1.8  meter  diameter,  3.7  meter  long  T6  vacuum  facility.  The  facility  is 
equipped  with  a  32  inch  diameter  diffusion  pump  with  an  estimated  pumping  speed  of  17,000  1/s  on  air.  The  facility 
is  also  equipped  with  cryogenic  pumps  with  an  estimated  pumping  speed  of  90,000  1/s  on  Xe. 

Facility  background  pressure  was  measured  with  an  INFICON  Bayard-Alpert  Pirani  Combination  Gauge, 
model  BPG400,  calibrated  for  xenon.  This  gauge  was  located  at  a  flange  on  the  side  of  the  test  facility 
approximately  2/3  of  a  meter  downstream  from  the  thruster  exit.  The  indicated  pressure  was  significantly  lower 
when  the  thruster  was  running  on  iodine  than  when  it  was  running  on  xenon  at  nearly  identical  conditions. 
However,  the  gauge  responds  differently  to  Xe,  I,  and  I2  and  the  correct  mix  is  not  known 

For  an  individual  gas,  the  actual  pressure  is  the  product  of  the  indicated  pressure,  pind  ,  and  a  calibration  factor, 
Cg .  The  BPG400  manual  identifies  calibration  factors  for  Xe  and  Kr,  but  not  for  I  or  I2.  The  relative  sensitivity, 
Rg ,  is  the  inverse  of  the  calibration  factor.  If  a  hot  ionization  gauge  has  not  been  calibrated  for  a  particular  gas,  the 
relative  sensitivity  may  be  estimated  from  electron  impact  ionization  cross  sections  at  150  eV,  e.g.  Rg  ~  o\,  /  <jN2 

Table  3  identifies  ionization  cross  sections  for  N2,  Kr,  Xe,  I  and  I2  from  the  literature”  ’  ’  ’  and  the 
corresponding  relative  sensitivities  with  respect  to  N2  and  Xe.  The  estimated  relative  sensitivities  for  I  and  I2  with 
respect  to  Xe  are  1.08  and  1.94,  respectively,  corresponding  to  calibration  factors  of  0.93  and  0.52.  Thus,  the 
indicated  pressure  (which  is  calibrated  for  Xe)  provides  a  very  conservative  estimate  of  the  actual  pressure  in  the 
facility.  The  estimated  relative  sensitivities  for  Kr  and  Xe  are  consistent  with  calibration  factors  in  the  BPG4000 
manual.  The  estimated  sensitivities  for  Kr,  Xe,  and  I2  are  also  consistent  with  empirical  sensitivity  factors  collected 
by  R.  Summers  at  NASA  Lewis  Research  Center.27 


Table  3  Calculating  the  Gas  Correction  Factor3 


Gas 

Ionization  Cross  Section 
at  150  eV  (Ang2) 

Calculated  Relative 
Sensitivity  w.r.t.  Xe 

Calculated  Relative 

Sensitivity  w.r.t.  N2 

Reported  Relative 
Sensitivity  w.r.t.  N2 

n2 

1.75 

0.40 

1.00 

Kr 

3.28 

0.75 

1.87 

<J) 

i 

Xe 

4.38 

1.00 

2.50 

2.2  -2.9 

1 

4.74 

1.08 

2.71 

h 

8.48 

1.94 

4.85 

5.4 

E.  Diagnostics 

Plume  density  was  measured  with  a  nude  Faraday  probe  originally  developed  at  MIT  which  has  been  cross- 
calibrated  against  a  larger  and  well-characterized  JPL  probe.28  29  The  Faraday  probe  was  mounted  on  an  automated 
rotating  arm  positioned  such  that  the  probe  swept  a  1 80  degree  arc  about  the  center  of  the  exit  plane  of  the  thruster, 
bisecting  the  plume.  The  pivot  axis  was  located  below  the  channel  exit,  and  the  radial  distance  from  the  probe  tip  to 
the  thruster  exit  was  approximately  37  cm. 

The  ExB,  ESA,  and  combined  probes  were  sourced  from  Plasma  Controls  LLC.30  These  probes  were  also 
mounted  upon  the  rotating  stage.  The  ExB  probe  makes  use  of  the  Lorentz  force  to  analyze  the  ion  population  by 
species.  However,  this  measurement  includes  ions  of  all  energies.  The  ESA  probe  uses  electrostatic  forces  to 
analyze  the  ion  population  by  charge  to  mass  ratio,  but  it  does  not  differentiate  between  species.  The  combined 
ESA/ExB  probe  is  able  to  analyze  the  ion  population  by  species  at  a  selected  energy  band.  The  combination  of 
measurements  provides  a  detailed  view  of  the  ion  energy  distribution. 

Thrust  was  measured  using  a  calibrated  thrust  stand  of  the  inverted  pendulum  type.31  However,  some  internal 
degradation  of  the  anode  was  observed  during  testing,  creating  uncertainty  as  to  whether  performance  changes  were 
related  to  beam  composition  or  to  non-uniformities  in  the  gas  flow  distribution.  Thus,  thrust  data  are  not  reported. 


3  The  ionization  cross  section  for  I2  (8.48  Ang2)  is  a  theoretical  value,  and  the  cross  section  for  I  calculated  by  the 
same  authors  using  the  same  methods  is  4.17  Ang2,  which  is  significantly  smaller  than  the  reported  experimental 
value  of  4.74  Ang2. 
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III.  Measurements 


A  Faraday  probe  was  initially  used  to  scope  out  the  plume  current  profde.  Then  the  ExB  probe  was  used  to 
obtain  an  overview  of  the  xenon  and  iodine  species  present  across  the  whole  plume  at  multiple  operating  conditions. 
Next,  the  ESA  probe  was  used  to  obtain  ion  energy  distributions  at  selected  operating  conditions.  Finally,  the 
combined  ESA/ExB  probe  was  used  to  determine  current,  mass  and  number  distributions  ( /  ,  g ,  and  h  )  within  the 
high  energy  populations  of  the  overall  ion  energy  distribution. 

A.  Faraday  Probe  Data 

At  the  beginning  of  testing,  the  beam  current  profde  at  several  conditions  was  measured  with  xenon.  The 
probe  was  swept  at  r  =  37  cm,  with  the  radial  axis  centered  roughly  underneath  the  thruster  exit.  The  angular 
increment  was  5  degrees,  and  an  offset  of  approximately  2  degrees  was  corrected  during  analysis.  The  probe  area 
was  corrected  per  the  method  of  Brown  and  Gallimore.  "  The  indicated  background  pressure  was  2.6xl0"5  torr  when 
the  facility  was  evacuated  with  a  diffusion  pump.  The  indicated  pressure  dropped  to  3.6xl0"6  torr  with  the  addition 
of  cryogenic  pumps.  This  had  a  marked  effect  upon  the  measured  current  density  profde,  particularly  in  the  wings. 

In  Fig.  7  the  current  density  at  Vd  =250  V  and  Pd  =200  W  with  cryo-pumps  operational  (which  is  roughly 
representative  of  iodine)  is  normalized  by  the  value  at  the  center  of  the  beam  and  multiplied  by  the  weighting 
factors  for  current  and  thrust,  sin($)  and  sin(#)cos(60 .  The  plot  of  j (9)  sin($)  cos($)  shows  that  the  axial 
momentum  is  almost  entirely  within  45  degrees  of  the  centroid,  but  peaks  at  12  degrees  -  not  at  the  center.  Viewing 
this  as  a  histogram,  approximately  18%  of  the  thrust  is  centered  around  12  degrees,  16%  is  centered  around  17 
degrees,  and  so  on.  Most  of  the  momentum  is  found  from  10  degrees  outward  -  a  region  where  significant  iodine 
dimer  populations  are  sometimes  observed. 


Fig.  7  Current  density  profile  for  thruster  flowing  Xe  with  weighting  factors 
B.  ExB  Probe  Data 

The  ExB  probe  was  used  to  obtain  an  overview  of  the  species  present  in  the  plume  as  a  function  of  9 ,  Vd  ,  and 
/ d  .  Sample  ExB  data  for  iodine  at  Vd  =250  V  and  Pd  =150  W  (Id  =0.6  A)  are  plotted  in  Fig.  8  for  9  <  50°.  This 

condition  was  selected  to  show  for  its  relatively  large  dimer  population,  seen  here  at  approximately  15  volts.  From 
left  to  right,  the  signals  represent  I2+,  I+,  I2  and  I.  Xenon  ExB  data  at  the  same  conditions  are  plotted  in  Fig.  9. 
From  left  to  right,  the  signals  represent  Xe+,  Xe2+,  and  Xe3+. 
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Fig.  8  Iodine  ExB  plume  sweep  at  250V  150W  (j)ind  »  4.5xl0'6  torr) 


< 

c 

OJ 

1- 

3 

u 

■o 

01 


5.E-10 


4.E-10 


3.E-10 


2.E-10 


<u  l.E-10 

5 


0.E+00 


10  15  20  25  30 

Set  Voltage,  V 


35  40 


- Odeg 

- 15  deg 

—  30  deg 

- 45  deg 

- 60  deg 


Fig.  9  Xenon  ExB  plume  sweep  at  250V  150W  (p,„ d  «2xl0'5  torr) 

The  size  of  each  population  was  determined  by  the  area  under  the  curve.  For  each  population,  area  A  was  taken 
to  be  the  height  of  the  peak  multiplied  by  the  width  at  half  maximum.  The  base  from  which  the  height  was 
measured  was  determined  by  averaging  the  low  energy  portion  of  the  data  (where  no  populations  exist).  The 
parameter  /  is  the  fraction  of  the  total  collected  current  by  species.  This  was  calculated  as 

fk=^r~-  (15) 

2> 

k 

For  iodine,  the  four  assumed  species  were  I2+,  I+,  I  ,  and  I  .  For  xenon,  the  species  were  Xe+,  Xe2+,  and  Xe3+. 
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At  0  <  20°,  ExB  measurements  indicate  /  associated  with  I2+  varies  between  0%  and  6%.  The  largest 
fractions  were  measured  at  low  discharge  current  ( I  d  =0.6  A).  Much  larger  I2+  fractions  were  measured  at  0  >  30°. 
However,  ESA  measurements  described  in  the  next  section  showed  that  low  energy  ions  are  a  significant  fraction  of 
the  total  beginning  at  0>\5a  .  To  determine  the  potential  effect  of  dimers  upon  performance,  the  low  energy 
portion  of  the  distribution  must  be  filtered  out.  This  is  accomplished  with  the  combined  probe. 

C.  ESA  Probe 


The  ESA  probe  was  used  to  determine  the  ion  energy  distribution  as  a  function  of  6  for  each  thruster  operating 
point  (  Vd  ,  Pd  )  of  interest.  Fig.  10  shows  the  iodine  ESA  sweep  at  Vd  =250  V  and  Pd  =150  W  (Id=0.6  A)  .  Fig. 

1 1  shows  the  xenon  ESA  sweep  at  the  same  discharge  conditions,  but  higher  background  pressure.  For  each  such  set 
of  sweeps  at  a  given  operating  condition,  the  significant  ion  populations  were  identified.  Near  the  beam  centroid, 
most  of  the  ions  were  highly  energetic.  Lower  energy  populations  increased  with  6  in  relative  magnitude.  At  0  > 
15,  the  velocity  distribution  is  far  from  Maxwellian. 

D.  Combined  Probe 


The  combined  ESA/ExB  probe  was  used  to  analyze  specific  populations  selected  from  the  ESA  data 
current  fractions  fk  were  determined  from  areas  as  before.  Number  fractions,  hk  ,  were  determined  by  di 
the  collected  current  associated  with  each  species  by  the  charge  associated  with  that  species: 

Ak  1 0  k 


K  = 


^Ak  /  ‘ Ik 


(16) 


k 

Mass  fractions,  gk ,  were  determined  by  multiplying  the  current  associated  with  each  species  by  mass,  and  dividing 
by  charge: 


o  _  /  qk 

^  AkM k  /  qk 

k 


(17) 


For  xenon,  h  and  g  are  equivalent. 


Iodine  ESA/ExB  data  for  Vd  =250  V  and  Pd  =150  W  ( Id  =0.6  A)  are  shown  in  Fig.  12.  Sweeps  were  taken  at 

the  primary  energy  peak  at  217V  for  all  angles.  The  indicated  tank  pressure  was  6.6x1 0'6  torr,  and  the  actual 
pressure  was  lower  by  a  factor  of  0.52  or  greater.  By  current,  the  measured  high  energy  dimer  fraction  at  the  beam 
centroid  was  2.2%  ( f k  «  0.022 ).  By  number  and  mass,  the  fractions  at  6  =0  were  2.3%  ( hk  ~  0.023  )  and  4.5%,  ( 
gk  =0.045).  Much  higher  fractions  were  measured  at  larger  angles.  For  instance,  at  0  =30  degrees,  the  dimer 
fraction  by  mass  was  22%  ( gk  >  0.20). 

Xenon  data  for  Vd=250  V  and  Pd  =150  W  ( I d  =0.6  A)  are  shown  in  Fig.  13  and  summarized  in  Table  5. 
Sweeps  at  250  V,  150  W  were  taken  at  primary  energy  peaks  of  215V  or  220  V.  The  indicated  background  pressure 
was  2.4xl0"torr.  By  current,  the  high  energy  Xe21  fraction  was  22%  ( f k  =0.22)  at  (9=30  degrees.  By  number 
and  mass,  Xe2+  was  13%  of  the  high  energy  plume  ( hk  =  gk  =  0.13)  at  this  same  location. 
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Fig.  12  Iodine  combined  ESA/ExB  plume  sweeps  taken  at  250  V,  150  W  (pinct  ~  6.6xl0'6  torr) 


Table  4.  Iodine  ESA/ExB  distributions  at  250  V,  150  W 


Angle  (deg) 

Species 

Mass 

Charge 

f 

g 

h 

0 

l2+ 

2 

i 

0.02 

0.04 

0.02 

0 

l  + 

1 

i 

0.92 

0.93 

0.95 

0 

.2+ 

1 

2 

0.05 

0.03 

0.03 

0 

.3+ 

1 

3 

0.00 

0.00 

0.00 

10 

i2+ 

2 

1 

0.04 

0.08 

0.04 

10 

r 

1 

1 

0.86 

0.87 

0.91 

10 

.2+ 

1 

2 

0.09 

0.05 

0.05 

10 

|3+ 

1 

3 

0.01 

0.00 

0.00 

20 

l2+ 

2 

1 

0.07 

0.14 

0.07 

20 

l  + 

1 

1 

0.80 

0.80 

0.86 

20  r 

1 

2 

0.13 

0.07 

0.07 

20 

.3+ 

1 

3 

0.00 

0.00 

0.00 

30 

l2+ 

2 

1 

0.12 

0.22 

0.12 

30 

r 

1 

1 

0.76 

0.72 

0.81 

30 

■2+ 

1 

2 

0.12 

0.06 

0.06 

30 

.3+ 

1 

3 

0.00 

0.00 

0.00 

40 

l2+ 

2 

1 

0.10 

0.18 

0.10 

40 

l  + 

1 

1 

0.85 

0.80 

0.87 

40 

|2+ 

1 

2 

0.05 

0.02 

0.03 

50 

l2+ 

2 

1 

0.11 

0.20 

0.11 

50 

l  + 

1 

1 

0.79 

0.75 

0.83 

50 

|2+ 

1 

2 

0.11 

0.05 

0.06 

60 

l2+ 

2 

1 

0.08 

0.17 

0.09 

60 

l  + 

1 

1 

0.80 

0.78 

0.85 

60 

|2+ 

1 

2 

0.12 

0.06 

0.06 

70 

l2+ 

2 

1 

0.00 

0.00 

0.00 

70 

1 

1 

0.91 

0.95 

0.95 

70 

|2+ 

1 

2 

0.09 

0.05 

0.05 
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Fig.  13  Xenon  combined  ESA/ExB  plume  sweeps  taken  at  250  V,  150  W 


Table  5  Xe  ESA/ExB  Data  at  250  V,  150  W 


Angle  (deg) 

Species 

Mass 

Charge 

f 

h=g 

0 

Xe+ 

1 

1 

0.89 

0.95 

0 

Xe2+ 

1 

2 

r 

0.08 

0.04 

0 

Xe3+ 

1 

3 

f 

0.03 

0.01 

15 

Xe+ 

1 

1 

F 

0.79 

0.89 

15 

Xe2+ 

1 

2 

r 

0.17 

0.10 

15 

Xe3+ 

1 

3 

r 

0.04 

0.01 

30 

Xe+ 

1 

1 

F 

0.75 

0.86 

30 

Xe2+ 

1 

2 

F 

0.22 

0.13 

30 

Xe3+ 

1 

3 

f 

0.03 

0.01 

45 

Xe+ 

1 

1 

F 

0.72 

0.85 

45 

Xe2+ 

1 

2 

f 

0.24 

0.14 

45 

Xe3+ 

1 

3 

f 

0.04 

0.01 

60 

Xe+ 

1 

1 

F 

0.71 

0.84 

60 

Xe2+ 

1 

2 

,F 

0.26 

0.15 

60 

Xe3+ 

1 

3 

F 

0.03 

0.01 

75 

Xe+ 

1 

1 

F 

0.82 

0.90 

75 

Xe2+ 

1 

2 

F 

0.18 

0.10 
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IV.  Discussion 


The  influence  of  multiple  species  upon  the  overall  efficiency  and  thrust  to  power  can  be  quantified  using  the 
molar  (number)  fraction  of  the  flux,  F.  Locally,  the  momentum  flux  is  proportional  to  ^hkMkvk  ,  where 

k 

vk  =  -y]2qkV /Mk  .  The  local  mass  flux  is  proportional  to  YjhkMk  •  If  1d  =IcxIb  ,  the  energy  cost  of  the  ions 

k 

at  constant  voltage  may  be  assumed  to  vary  in  proportion  to  the  ion  current  density  and,  this,  in  proportion  to 
Y \hkqk  .  (In  actuality,  with  each  of  the  k  species,  there  may  be  a  wide  variation  in  point  of  origin).  With  these 

k 

assumptions,  in  analogy  to  equation  (7),  a  species  efficiency,  i)s ,  may  be  written: 

(Ya^v,)2 

Vs  «  - • 

2_,hkMk  2jhVk 

k  k 

In  a  similar  fashion,  a  thrust  to  power  factor  may  be  written: 

kM  kvk 

hkclk 
k 

For  the  three  primary  ionic  species  in  the  beam,  normalized  values  of  charge,  mass,  velocity,  momentum,  and 
energy  are  shown  in  Table  6.  Although  I21  is  twice  as  energetic  as  land  I2+,  each  I2  ion  contributes  two  electrons  to 

the  discharge  instead  of  one.  I2  is  V2  times  as  fast  as  I+  and  twice  as  fast  as  I2+,  but  I2+  and  F  have  the  same 
momentum.  With  these  values,  relative  values  of  rjs  and  TIP  can  be  estimated  from  equations  13  and  14  as  /?, , 
h2  ,  and  h2  are  varied. 

Table  6  Iodine  Species  and  their  Relative  Properties 


(13) 


(14) 


Parameter 

Symbol 

r 

i2+ 

l2+ 

Charge 

q 

i 

2 

1 

Mass 

M 

i 

1 

2 

Charge/Mass 

q/M 

i 

2 

0.5 

Energy 

8 

i 

2 

1 

Velocity 

V 

i 

1.41 

0.71 

Momentum 

p 

i 

1.41 

1.41 

Momentum/Energy 

p/s 

i 

0.71 

1.41 

Fig.  14  plots  Equation  13  as  a  function  of  the  singly  charged  dimer  (I2+)  and  monomer  (I  ),  fractions,  assuming 
the  rest  of  the  discharge  is  comprised  of  doubly  charged  monomers  (I2+).  The  curves  represent  I+  fractions.  The  X- 
axis  is  the  l2  fraction.  The  presence  of  multiple  species  always  exacts  an  efficiency  loss.  However,  rjs  is  always 
greater  than  93%  when  more  than  50%  of  the  beam  is  I+.  This  result  also  is  significant  for  its  relation  to  two 
thrusters  operating  next  to  each  other  at  different  specific  impulse,  or  a  multiple  channel  “nested”  Hall  thruster 
where  the  two  channels  operate  at  specific  impulse;  a  modest  difference  in  specific  impulse  has  only  a  minor  impact 
upon  efficiency. 

Fig.  15  plots  the  net  thrust  to  power  (equation  14)  as  function  of  the  I2+  and  I+  fractions,  again  assuming  the  rest 
of  the  discharge  is  comprised  of  I.  This  figure  shows  that  I2+  could  increase  pis  significantly. 
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Fig.  14  Efficiency  Factor  as  a  Function  of  Molar  (Number)  Fractions 
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Fig.  15  Thrust  to  Power  Factor  as  a  Function  of  Molar  (Number)  Fractions 


The  following  example  helps  clarify  the  meaning  of  the  plots.  If  the  flux  is  50%  I+  ( =  0.5) ,  the  plots  assume 
h2  +  h2  =  0.5  .  The  worst  case  for  efficiency  occurs  when  h2  =  h2  (equal  molar  fractions  of  I2+  and  l2  ),  at  which 
i]s  =  0.933  (Fig.  14)  and  pi  s  =0.966  (Fig.  15).  The  best  case  for  efficiency  occurs  when  h2  =  0  or  h2  =  0 ,  in 
which  case  i]s  =  0.971.  In  this  best  case,  if  there  is  no  I2+,  then  p/s=  1.21.  If  there  is  no  I2+,  then  /?  /  £•  =0.805. 

The  ExB  probe  measurements  showed  the  presence  of  molecular  species  throughout  the  central,  momentum 
carrying  portion  of  the  plume.  However,  the  dimer  population  varied  with  operating  conditions  and  seemed  to  be 
concentrated  off  axis.  Table  7  summarizes  the  ESA/ExB  current  fractions,  /  ,  across  a  broad  range  of  operating 
conditions.  These  data  seem  to  indicate  that  the  dimer  fraction  is  largest  when  the  current  is  low  and  the  discharge  is 
rarified.  The  highest  dimer  fractions  by  current  are  observed  at  250  V  and  150  W,  where  1=0.6  A.  Relatively  large, 
localized  populations  were  also  found  at  Id  =0.71  and  I d  =0.8  A. 
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Table  7  Summary  of  iodine  F.SA/F.xB  Prohe  Current  Fractions  at  Hinh  F.nerpy  Peaks 


Discharge  Conditions 

150 

200 

250 

250 

300 

350 

350  [V] 

Angle  [deg] 

Species 

r  1.00 

r  1.00 

r  0.60 

r  0.80 

r  0.80 

r  0.71 

r  1.00  [A] 

150 

200 

150 

200 

240 

250 

350  [W] 

0 

lz+ 

0.02 

0.01 

0.01 

0 

f 

0.53 

1.00 

0.92 

0.95 

0.93 

0.94 

0.96 

0 

■  2+ 

0.47 

0.00 

0.05 

0.04 

0.07 

0.05 

0.04 

0 

.3+ 

0.00 

0.00 

0.00 

10 

i2+ 

0.04 

0.03 

0.02 

0.02 

10 

f 

0.98 

0.83 

0.86 

0.82 

0.83 

0.90 

0.88 

10 

■  2+ 

0.02 

0.14 

0.09 

0.08 

0.11 

0.08 

0.12 

10 

.3+ 

0.03 

0.01 

0.07 

0.05 

20 

i2+ 

0.07 

0.02 

0.02 

0.06 

0.00 

20 

f 

0.99 

0.40 

0.80 

0.85 

0.82 

0.77 

0.86 

20 

0.01 

0.60 

0.13 

0.13 

0.15 

0.18 

0.13 

20 

.3+ 

0.00 

0.01 

0.00 

30 

i2+ 

0.12 

0.00 

0.03 

0.02 

30 

f 

1.00 

0.85 

0.76 

0.95 

0.77 

0.86 

0.79 

30 

l2+ 

0.15 

0.12 

0.00 

0.18 

0.14 

0.19 

30 

.3+ 

0.00 

0.05 

0.03 

40 

i2+ 

0.10 

0.00 

0.05 

0.01 

0.02 

40 

f 

1.00 

0.85 

0.85 

0.90 

0.76 

0.87 

0.75 

40 

■  2+ 

0.15 

0.05 

0.10 

0.19 

0.13 

0.23 

50 

i2+ 

0.11 

0.00 

0.04 

0.01 

0.00 

50 

f 

1.00 

1.00 

0.79 

0.91 

0.82 

0.90 

0.84 

50 

■  2+ 

0.11 

0.09 

0.14 

0.09 

0.16 

60 

i2+ 

0.08 

0.06 

0.04 

0.08 

0.00 

60 

f 

1.00 

0.80 

0.81 

0.92 

0.81 

0.83 

60 

■  2+ 

0.12 

0.13 

0.04 

0.12 

0.17 

70 

i2+ 

0.00 

0.04 

70 

f 

1.00 

0.91 

0.95 

1.00 

0.81 

70 

i2+ 

0.09 

0.00 

0.19 

Fig.  16  plots  the  mass  fraction,  gk,  measured  verses  angle  for  the  250  V,  150  W  conditions.  The  high  energy 
dimer  fraction  is  smaller  near  the  beam  centroid  than  off  axis.  The  measured  i2+  fraction  is  greatest  from  30  to  50 
degrees,  which  is  the  peak  of  the  weighting  function  for  thrust  given  in  equation  1 1 .  Dimers  are  also  present  in  the 
far  field  (50-60  degrees),  but  these  diminish  in  significance  due  to  the  exponential  decrease  in  current  density 
compounded  by  the  as  the  weighting  factor  cos(Y/)  sin( 0) . 

The  measurements  indicate  that  the  dimer  fraction  at  some  conditions  is  sufficient  to  increase  thrust  to  power  by 
a  small,  but  measurable  amount.  Per  Table  4,  at  250  V,  150  W  the  data  show  a  dimer  number  fraction  between  2 
and  12  percent  in  the  central  part  of  the  plume.  Therefore,  per  Fig.  14  the  species  or  poydispersivity  efficiency  is 
likely  to  be  98  to  99%  and  the  thrust  to  power  factor,  pis  will  be  between  0.95  and  1.05,  depending  on  the  detailed 
distribution  of  momentum  vs.  angle. 
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Fig.  16  Mass  fraction  of  flux  at  250  V,  150  W 


This  method  used  to  determine  /  is  similar  to  the  “method  of  triangle  fitting”  recommended  by  Shastry  et  al. 

and  suggested  by  Beal."  ’1 2 3 4 5 6  However,  Shastry  recommends  multiplying  the  peak  height  by  the  half  width  at  half 
maximum,  where  the  half  width  is  defined  on  the  high  energy  side  of  the  curve.  The  analysis  here  uses  the  full 
width.  To  determine  the  sensitivity  of  the  results  to  this  difference,  a  subset  of  the  iodine  data  at  250  V,  150  W  were 
re-analyzed  using  both  methods.  In  relative  terms,  using  the  half  width  on  the  high  side  decreased  the  dimer  current 
fractions  by  9,  12,  and  20%  at  20,  30,  and  40  degrees,  respectively.  In  absolute  terms,  the  difference  was  0.5,  1.4 
and  1.7%.  Shastry  also  recommends  the  use  of  a  charge  exchange  correction  model.  No  such  correction  was 
applied  to  the  data  because  the  composition  of  the  background  gas  is  unclear. 

V.  Conclusions 

Prior  performance  measurements  showed  that  iodine  could  be  an  attractive  replacement  for  xenon  in  Hall  effect 
thrusters.  Overall  performance  was  found  to  be  on  par  with  xenon,  while  initial  plume  measurements  showed  a  well 
collimated  beam  profile  and  a  small  but  significant  population  of  ionized  dimers.  At  some  conditions,  higher  thrust 
to  power  was  measured  with  iodine  than  with  xenon. 

The  data  reported  here  lend  further  insight  into  the  plume  structure  and  ion  species  distribution.  Significant 
populations  of  diatomic  iodine  were  observed  in  the  plume.  These  dimers  were  concentrated  slightly  off  axis  from 
the  thruster  centroid,  and  were  observed  to  vary  with  thruster  operating  conditions.  Such  dimers  could  change  the 
observed  thruster  efficiency  and  thrust  to  power  by  measurable  amounts.  In  any  case,  dimers  should  be  accounted 
for  in  future  analytical  and  numerical  studies  of  iodine  fueled  thrusters. 
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